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We investigate the entanglement properties of the Kondo spin chain when it is prepared in its ground 
state as well as its dynamics following a single bond quench. We show that a true measure of 
entanglement such as negativity enables to characterize the unique features of the gapless Kondo 
regime. We determine the spatial extent of the Kondo screening cloud and propose an ansatz for 
the ground state in the Kondo regime accessible to this spin chain; we also demonstrate that the 
impurity spin is indeed maximally entangled with the Kondo cloud. We exploit these features of 
the entanglement in the gapless Kondo regime to show that a single local quench at one end of a 
Kondo spin chain may always induce a fast and long lived oscillatory dynamics, which establishes 
a high quality entanglement between the individual spins at the opposite ends of the chain. This 
entanglement is a footprint of the presence of the Kondo cloud and may be engineered so as to 
attain - even for very large chains- a constant high value independent of the length; in addition, it 
is thermally robust. To better evidence the remarkable peculiarities of the Kondo regime, we carry 
a parallel analysis of the entanglement properties of the Kondo spin chain model in the gapped 
dimerised regime where these remarkable features are absent. 

1 Introduction 

Kondo systems ||2] [H are expected to be very distinctive in the context of entanglement for at least two 
reasons: (a) Despite being "gapless", they support the emergence of a length scale i§ EKll which should 
be reflected in the entanglement, making it markedly different from that in the more conventional gapless 
models studied so far; (b) They are expected to have a more exotic /orm of entanglement than the widely 
studied spin-spin and complementary block entanglements. Indeed, in Kondo systems, the impurity spin 
is expected to be mostly entangled with only a specific block of the whole system. This is, of course, 
merely an intuition which needs to be quantitatively verified with a genuine measure of entanglement: 
this task has been accomplished in |f3l where we provided the only characterization of the Kondo regime 
based entirely on a true measure of entanglement such as negativity H. 

The simplest Kondo model El |5l describes a single impurity spin interacting with the conduction 
electrons in a metal; the ground state is a highly nontrivial many body state in which the impurity spin 
is screened by conduction electrons in a large orbital of size - the so called Kondo screening length. 
Many physical observables vary on the characteristic length scale ^ , which is a well defined function of 
the Kondo coupling 121 . The screening length ^ determines the spatial extent of the Kondo cloud whose 
signatures in physical systems have been so far a challenging problem repeatedly addressed by various 
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means |[Tl|6l|7l. As we shall detail in the following an entanglement measure enables to fully determine 
i§ 111 and to use its knowledge to engineer long range distance independent entanglement HI. 

Recently [9], it has been pointed out that the universal low energy long distance behavior of this 
simple Kondo model arises also in a spin chain when a magnetic impurity is coupled to the end of a 
gapless Heisenberg anti-ferromagnetic Ji — Ji spin 1/2 chain. The spin chain Kondo model ||9l is defined 
by the Hamiltonian 

Hj = J'{Jia\.02+J20i.Oj) 

N-\ N-2 
i=2 i=2 

where a,- = {of, of, of) is a vector of Pauli operators at site /, N is the total length of the chain, J2 is the 
next nearest neighbor coupling and the nearest neighbor coupling J\ is normalized to 1. The impurity 
spin, located at one end of the chain, is accounted for by weaker couplings to the rest of the system; in 
the following, see Fig. [T^, both couplings J\ and J2 are weakened by the same factor J' , which quantifies 
then the impurity strength, for < ^2 < -/f = 0.2412, the spin system is gapless and it supports a Kondo 
regime ifTOl [Tl. For J2 > /f' the system enters the gapped dimer regime, where the ground state takes 
a dimerised form; at the Majumdar-Ghosh fm point {J2 = 0.5), the ground state becomes just a tensor 
product of singlets. For J2 > 0.5, incommensurability effects |[T2| emerge. 

In the following we detail our analysis ||3l [H of the remarkable features of entanglement in the 
Kondo regime 111 [H of the spin chain Kondo model. First of all, we use negativity Q to characterize 
|[3l the entanglement in the Kondo regime; namely,for this spin chain in the Kondo regime: (i) we 
demonstrate that the impurity spin is maximally entangled with the Kondo cloud; (ii) we determine the 
spatial extent of the Kondo screening length using only an entanglement measure; (iii) we motivate an 
ansatz for the ground state in the Kondo regime; (iv) we evidence the scaling of negativity as pertinent 
parameters are varied. Then, we use the knowledge acquired in 131 to engineer- for chains of arbitrary size 
N- a long range distance-independent entanglement through a non-perturbative dynamical mechanism 
which requires only a minimal action on a spin chain, namely a sudden quench of a single bond |[8l . 
The ensuing long range distance independent entanglement- as well as the remai^kable entanglement 
oscillations observed in our numerical simulations [8]- are a footprint of the emergence of the length 
scale ^ in an -otherwise gapless- system. To accomplish these tasks we designed fJl a Density Matrix 
Renormalization Group (DMRG) approach enabling to investigate the entanglement between a single 
spin and a pertinent block of the chain. We also developed HI a time dependent DMRG to simulate the 
dynamics of the system after a sudden quench of a local bond; to better evidence the unique properties 
of the entanglement in the Kondo regime we carried a parallel analysis of the entanglement properties 
of this model in the gapped dimerised regime where all these remarkable features of entanglement are 
absent. 

The investigation of entanglement in many-body condensed matter systems is cun^ently a topic of 
intense activity |[T3l[T4l[T5l[T6l[r7l[T8l[T0l [Tl[T9l[3l[8l. In many instances investigations have focussed 
on the entanglement between individual elements, such as single spins, or the entanglement between 
two complementary blocks in the ground state of a condensed matter system. The former is generically 
non-zero only between nearest or next to nearest neighbors |[T4l IT5l ; at variance, for complementary 
blocks, the whole system is in a pure state and the von Neumann entropy is a permissible measure of 
the entanglement. In conventional gapless phases, due to the absence of an intrinsic length scale, the 
von Neumann entropy diverges with the size of the blocks |[T6l [TTl [TSl ; at variance, our analysis |[3l [H 
provides for the first time a characterization of the entanglement in gapless regimes of a many body 
system where a length scale emerges as a result of the presence of a Kondo impurity. 
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From a more practical viewpoint there is now a huge demand for long range entanglement between 
individual spins in quantum information theory. For instance, the long range entanglement between in- 
dividual particles gives the opportunity of implementing teleportation [1201 which offers perfect quantum 
communication between distant parties. Though entanglement in condensed matter systems is typically 
very short ranged llT5l . there have been a few other proposals for long distance entanglement; however, 
they come at a high price. For instance, there are proposals exploiting weak couplings of two distant 
spins to a spin chain Il2ni22l . which have very limited thermal stability or very long time-scale of entan- 
glement generation. Otherwise, a dynamics has to be induced by lai^ge-scale changes to the Hamiltonian 
of a system |[23l . 

The paper is organized as follows: In section ^ we review our analysis 131 of the static entanglement 
properties of the ground state of the Kondo spin chain model using negativity lH as a true measure of 
entanglement: there, using only quantum information tools, we are ableQ to determine the spatial extent 
of the Kondo cloud as well as to provide an ansatz for the ground state in the Kondo regime of the spin 
chain Kondo model. In section ^ we point out fSl how - only in the Kondo regime of the spin chain 
Kondo model- a local quench of the last bond of the Kondo spin chain induces a thermally robust high 
quality long range and distance-independent entanglement between the two individual ending spins of 
the chain: there, we emphasize how the onset of this long range distance-independent entanglement may 
be regarded as the footprint of the emergence of the Kondo cloud in the Kondo regime of the spin chain 
Kondo model. Section ^ is devoted to a short discussion of our results. 

2 Entanglement in the Ground State of the Spin Chain Kondo Model 

In this section we detail the results of our analysis IJl of the entanglement properties of the ground state 
of the Hamiltonian described by ([T])- We show that a true measure of entanglement enables one to use 
only quantum information tools to determine the Kondo screening length ^ and to provide an ansatz for 
the ground state of this spin chain in the Kondo regime. 

A true measure of entanglement should satisfy a set of postulates - for example, it should be non- 
increasing under local actions: such a genuine measure does exist for two sub-systems of arbitrary size 
even when their combined state is mixed, as it happens in Kondo systems. This measure is the negativity 
im and it has been successfully used to quantify the entanglement in a harmonic chain |[24l l25l and 
between distant regions of critical systems |[26ll27l . For bipartite systems, negativity is 

defined as £ = — 1, where a, denote the eigenvalues of the partial transpose of the whole 

density matrix of the system with respect to one of the two subsets of the given partition and | . . . | is the 
absolute value H. 

In the following we use DMRG as the numerical tool needed to measure the static entanglement 
properties observed in the Kondo regime of the spin chain Kondo model. 

2.1 The Entanglement Healing Length 

To study the entanglement of the ground state we divide- see Fig. \Vp - all the spins of the chain in 
three different groups: the impurity spin, block A, which contains the L spins next to the impurity (L = 
0, 1, ...,A'^ — 1) and block B formed by the remaining N — L — I spins. We used |[3l negativity to fully 
characterize the entanglement between the impurity and block B in both the gapless Kondo and the 
gapped dimerised regimes. 

We determined IS the size of the block A when the entanglement between the impurity and block 
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Figure 1 : (a) Kondo Spin chain with next nearest neighbor Heisenberg interaction with one impurity at 
one end. (b) The chain is divided into three parts, an impurity, a block A and a block B. Entanglement is 
computed between the impurity and block B. 



B is almost zero; by this procedure one defines |3] an Entanglement Healing Length (EHL) L*, i.e. the 
length of the block A which is maximally entangled with the impurity. It should be emphasized that 
an EHL may be defined for all spin chains: however, we showed |[3l that, only in the gapless Kondo 
regime of the Kondo spin chain, EHL scales with the strength of the impurity coupling just as the Kondo 
screening length, ^ , does. 

In the gapless regime of the Kondo spin chain, measuring the EHL through negativity [3] yields a 
genuine quantum information tool to detect the Kondo screening length |[TJ|6l|7l. In addition v/e found 
that- in the Kondo regime- entanglement, as quantified by negativity, is a homogeneous function of the 
two ratios: N/L* and L/N, where L is the size of the block A, i.e. the block adjacent to the impurity, 
and A'^ is the length of the whole chain. As a result, the entanglement in the Kondo regime is essentially 
unchanged if one rescales all the length scales with the EHL L* . Of course, EHL can be defined also in 
the gapped dimerised regime but negativity is now a function of the three independent quantities A'^, L 
and L*: i. e., scaling of entanglement with EHL is absent in the gapped dimerised regime 13J. 

2.2 Numerical Approach: Density Matrix Renormalization Group 

The DMRG ||28| approach has been used Q to compute the ground state of the spin chain Kondo model 
for large chains up to A'^ = 250; in order to avoid finite size effects, we took A'^ to be even and, thus, 
avoided problems arising from accidental degeneracies. 

In a DMRG approach the ground state of the system is partitioned in terms of states of a left block, 
a right block (not to be confused with blocks A and B) and two intermediate spins as shown in Fig. [2^. 
The states of the intermediate spins ai^e given in the computational (| t), | \.)) basis, while the states of 
the both blocks are usually in some non-trivial truncated DMRG basis. In this approach one has several 
representations for the ground state which vary due to the number of spins in the left (right) block and it 
is possible to go from one representation to the other by applying pertinent operators on each block. 

The main issue in the DMRG is that the dimension of the left (right) block is kept constant and in- 
dependent of whatever spins are there in that block. To have a fixed dimension for the left (right) block 
we truncated the Hilbert space so that the amount of entanglement between the two parts of the chain 
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Figure 2: (a) DMRG representation of the state of the chain keeps two intermediate spins in ordinary 
computational basis and the left and right blocks in a truncated DMRG basis, (b) The intermediate spin 
next to the impurity is traced out from the density matrix of the chain. This tracing is equivalent to adding 
the traced out spin to block A. (c) The basis of the right block of DMRG representation is transformed 
so that a single spin in the left side of the right block is represented in the computational basis while the 
state of the new right block is given in a DMRG basis. 



remains almost unchanged ||28l . To improve the precision of our results we swept |[3l all representations 
of the ground state for several times to get the proper basis for the left and the right blocks of all represen- 
tations. After some sweeps, when the ground state energy converges (we kept states for which the eiTor 
on the energy is less than 10^^), we paused to compute the entanglement. We took a representation of 
the ground state in which the left block contains just the single impurity spin and the right block contains 
A'^ — 3 spins: as a result, the single impurity spin is given in the computational basis and this allowed us 
to compute the negativity later. 

From the DMRG state, one should trace out the spins belonging to block A before computing the 
entanglement between the impurity and block B since it is most convenient to compute the entanglement 
between the impurity and the block B: due to the entanglement monogamy, this provides an equivalent 
information about the entanglement of the impurity with the block A. Our tracing procedure started with 
the density matrix of the ground state of the system in the representation shown in Fig. |2^; at this stage, 
the number of spins in the block A is zero (no spin has been traced out), all spins except the impurity 
belong to the block B, and the entanglement between the impurity and the block B is maximal (i.e. 
E = 1). Then, we traced out the intermediate spin next to the impurity as shown in Fig. |2Jd; this amounts 
to putting that spin into block A. Finally, as shown in Fig. |2j:, we transformed the DMRG basis of the 
right block so as to put a single spin at the left of the right block in the computational basis, while the 
state of the new right block is given in a DMRG basis. As a consequence, the resulting density matrix 
had the exact form of Fig. |2h and we continued the procedure to trace one spin at each step (i.e., put 
more spins in the block A) and computed the entanglement between the impurity and block B. 

2.3 Scaling of Negativity and Ansatz for the Ground State in the Kondo Regime 

As evidenced in |[3l, there is an EHL L* so that, for L> L* , the entanglement between the impurity and 
block B is almost zero: L* provides us with an estimate of the distance for which the impurity is mostly 
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Figure 3: (a) L* vs. for both Kondo (Jj = 0) and dimer regime {J2 = 0.42). (b) Entanglement vs. 

L/N for fixed N/L* = 4 when Jj = 0. (c) Entanglement vs. L/N for fixed N/L* = 4 at the critical point 
h = Ji- (d) Entanglement vs. L/N for fixed N/L* = 4 in the dimer regime {J2 = 0.42). 



entangled with the spins contained in block A. For large chains (A'^ > 200) in the Kondo regime, one finds 
that L* is almost independent of A'^ and depends only on /. In the Kondo regime, i.e. for J2 < J2, L* 
depends on J' just as the Kondo screening length ^ does HOlSl; for small /, L* oc e''/^, where a IS a 
constant. We plot L* as a function of 1/ VJ' in Fig. [3^. In a semi-logarithmic scale, the straight line plot 
exhibited in the Kondo regime {J2 = 0) shows that L* may be indeed regarded as the Kondo screening 
length. Moreover, the nonlinearity of the same plot in the dimer regime {J2 = 0.42), especially for small 
J', shows that, in the gapped dimerised regime accessible to the Kondo spin chain model, no exponential 
dependence on 1 / y/f holds. 

There is 131 also a remarkable scaling of negativity in the Kondo regime. This scaling may be 
regai^ded as yet another independent evidence of the fact that L* is indeed the Kondo length ^ . In general, 
the entanglement E between the impurity and block B is a function of the three independent variables, 
J',L and which, due to the one to one correspondence between J' and L*, can be written as E{L* ,L,N). 
We find that, in the Kondo regime, E = E{N /L* ,L/N). To illustrate this, we fix the ratio N/L* and plot 
the entanglement in terms of L/N for different values of /' (or equivalently L*) for /i = (Fig. [3J5) and 
for J2 = J2 (Fig- Et)- The complete coincidence of the two plots in Figs. [3J5 and c shows that, in the 
Kondo regime, the spin chain can be scaled in size without essentially affecting the entanglement as long 
as L* is also scaled. In the dimer regime the entanglement stays a function of three independent variables, 
i.e. E = E{L* ,L,N), and, as shown in Fig. [3}l, the entanglement does not scale with L*. In our approach, 
the EHL L* may be evaluated in both the Kondo and the dimer regime: the scaling behavior, as well 
as the dependence of L* on J', discriminates then between the very different entanglement properties 
exhibited by the spin chain Kondo model as J2 crosses Jj. 

We defined L* such that there is no entanglement between the impurity and block B when block A is 
made of L* spins. Conventional wisdom based on previous renormalization group analysis suggests that, 
in both regimes, the impurity and the block A of length L* form a pure entangled state, while block B 
is also in a pure state. This is indeed approximately true in the dimer regime (exactly true for J2 = 0.5) 
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Figure 4: (a) Entanglement vs. L for the two fixed points / — > and /' — 1 in tiie Kondo regime J2 = 0. 
(b) Exponential decay of entanglement in terms of L in a chain of length N = 250 for J' = 0.6. (c) 
Non-exponential decay for small /' in the dimer regime J2 = 0.42. 



but it turns out to be dramatically different in the Kondo regime. To check this, we computed the von 
Neumann entropy of the block B when block A has L* spins and found it to be non zero. Thus, the 
blocks A and B are necessarily entangled in the Kondo regime as there is no entanglement between the 
impurity and B. In fact, after a distance L*, the impurity is "screened" i.e, the block B feels as if it is part 
of a conventional gapless chain and has a diverging von Neumann entropy. The Kondo cloud is, then, 
maximally entangled with the impurity as well as being significantly entangled with block B. Based 
on the above, a simple ansatz for the ground state \GS) in the Kondo regime has been proposed in lO; 
namely, one may conjecture that 

where are constants, {|Lt (/')), \Lj{J'))} and are sets of orthogonal states on the cloud and 

the remaining system, respectively. At the fixed point /' — )• all spins except the impurity are included 
in and \Lj{J')). At / 1, very few spins are contained in and \Lj{J')) while {\Ri{J'))} 

represents most of the chain. 

For what concerns the mere evaluation of the amount of entanglement as J' is varied, we plot, in 
Fig. Ilh, the negativity as a function of L near by the two fixed points, i.e. 7' — and /' — 1, accessible 
in the Kondo regime: as expected, near /' — (i.e, for large values of the Kondo screening length), 
the entanglement remains large for rather large values of L while, for / — 1 (i.e. for a very small 
cloud) it decreases rapidly with L. Fig. |4^ (semi-logarithmic) shows that, also at the extreme limits 
7' — )• and 7' — )• 1, the entanglement decays exponentially with L since this a characteristic mark of 
the entanglement in the Kondo regime. This exponential decay of entanglement is absent in the dimer 
regime: Fig. |4j) shows that, in the gapped dimerised regime, only for rather large the entanglement 
decays exponentially with L while, for small /', the entanglement between the impurity and block B 
decays slower than an exponential as a function of L exhibiting even a plateau at short distances. This 
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latter feature is evidenced in Fig. |4]:, and is consistent with the emergence, for small J', of long range 
valence bonds between the impurity and far spins as a consequence of the onset of the dimerised ground 
state im. In fact, when J' is small, J2J' becomes much less than and the impurity forms valence 
bonds with distant spins, while the other spins, since for them J2 > /f' form a valence bond with their 
nearest neighbor to preserve the dimerised nature of the ground state: as a result, the impurity shares less 
entanglement with nearby spins and fulfills its capacity of entanglement forming valence bonds with the 
more distant spins in the chain. 



3 LONG-RANGE DISTANCE-INDEPENDENT ENTANGLEMENT IN 
THE KONDO REGIME 



Entanglement between a single impurity spin and a group of spins- such as the ones inside the Kondo 
cloud emerging in the Kondo regime accessible to the spin chain Kondo model- cannot be used as a 
resource for computational tasks since manipulation of many particles is an extremely difficult task. It 
is much more convenient, instead, to use entanglement between distant individual particles since, due to 
their localization, one may more easily resort to unitary gates and measurements to control these indi- 
vidual spins. In lU we proposed a procedure to convert the useless entanglement between the impurity 
and the cloud into the useful entanglement between the two ending spins of the Kondo spin chain. 

To engineer the long range entanglement between distant individual spins we took lH the finite 
Kondo spin chain ([T]) in its ground state \GSi) and, then, pertinently quenched the coupling at the op- 
posite end of the impurity allowing for the dynamics to develop entanglement. We showed that, in the 
Kondo regime, the entanglement between the two ending spins oscillates between high peaks with a 
periodicity determined by J' , while the dynamics is very fast (thereby decoherence hardly gets time to 
act) and is robust against thermal fluctuations. At variance, in the gapped dimer regime, the dynamics 
is much slower, qualitatively different and, in finite chains, it still generates some entanglement due to 
the unavoidable end to end effects, which are drastically tamed if one "cuts off the impurity from the 
chain during the dynamics. In the Kondo regime, cutting off the impurity has minimal effect on the 
final entanglement between the ending spins since, here, the process is induced by the emergence of the 
Kondo cloud H. 

3.1 Entanglement Oscillations Induced by Local Quench Dynamics in the Kondo Regime 

Initially, the system is assumed to be in the ground state \GSi) of Hj. A minimal quench modifies only 
the couplings of the A'^th spin by the amount /' (same as /' in Eq. ([T])) so that Hj is changed to 



Since \GSi) is not an eigenstate of Hp it evolves as \^{t)) = e^'^'''\GSj) . An entanglement E{N,t,J') 
between the ending spins emerges as a result of the above evolution. 

To compute E{N,t,J'), in [H we computed the reduced density matrix piN{t) = tfiy\¥{f)){¥{t)\ of 
spins 1 and N by tracing out the remaining spins from the state \ Y{t))- Then, we evaluated E{N,t,J') in 
terms of a measure of entanglement valid for arbitrary mixed states of two qubits called the concurrence 
||29l . We showed that entanglement took its maximum value E,„ at an optimal time tgpt and at optimal 



Hp 




(3) 
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Figure 5: Comparing the Kondo {J2 = 0) and dimer {J2 = 0.42) regimes, (a) Entanglement vs. time 
for N = 30 with /' = 0.19 for the Kondo regime (solid line) and J' = 0.44 for the dimer regime (dashed 
line), (b) £,„ vs. the length A^. (c) topt vs. A^. (d) 7',^^ vs. length A^. 



coupling J'^jpi such that Em = E{N,topt,J'opt)- We evidenced in HI that 7^^, is not a perturbation of /i and 
J2- If, as expected from scaling in the Kondo regime |[30l l3l. the dependence on A'^ and t can be accounted 
for by a redefinition of J' (equivalently <^), then tgpt and 7^pr cannot be independent quantities. 
Our numerical analysis showed that, in the Kondo regime, tQpf oc 

N and that J'^p; yields ^ =N-2; 

since ^ oc one gets 

t^,pf oc A^ oc e«/ V^' 

. For our choice of J',Ji and J2 the spin-chain dynamics is not analytically solvable and one has to resort 
to numerical simulations. Here, for A'^ > 20, we used the time-step targeting method, based on the DMRG 
algorithm introduced in ||3T1 . For A'^ < 20, instead, one may resort to exact diagonalization. 

We focused only on the first period of the entanglement evolution in both regimes, since both deco- 
herence and numerical errors make it unwise to wait for longer times. Fig. [Sja) shows that fast long-lived 
(non-decaying) periodic oscillations with a period of Itopt characterize the time evolution in the Kondo 
regime and that the maximal entanglement is achieved when the impurity coupling / equals the value 
J'gpf associated to a Kondo cloud of size ^ = N — 2 (the Kondo cloud generated by the impurity sitting 
on the left side touches the other side of the chain); in the dimer regime the dynamics appears more dis- 
persive and not oscillatory for any J'. In Fig. [^b) we plot the maximum of entanglement, E,„, induced 
by bond quenching as a function of the length A'^: though the entanglement decreases as A'^ increases, its 
value, in the Kondo regime, stays rather high and becomes almost distance independent for very long 
chains; furthermore, as A'^ increases, the entanglement generated in the Kondo regime is significantly 
bigger than the one in the dimer regime. 

Despite its lower value, achieving entanglement in the dimer regime costs more time, as shown in 
Fig- Etc). It is also cleai^ from Fig. [Sjc) that topt increases by A'^ linearly. At variance, in the Kondo 
regime, J'^p, slowly decreases as A'^ increases while it stays essentially constant in the dimer regime (Fig. 

md)). This is commensurate with the expectation that topi °^ e^^V^ in the Kondo regime, while, in the 
dimer regime, topt and 7^^, are two independent quantities. The plot in Fig. [6^ shows the exponential 

/— 1 /2 

dependence of tgpt on Jopt realized, in the Kondo regime, for long enough chains. 
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3.2 Probing the Size of the Kondo Cloud through Entanglement's Quench Dynamics 

How the dynamics allows to engineer - even for very long chains of size A'^- high entanglement oscilla- 
tions between the ending spins of the chain in the Kondo regime? To understand this, one should recall 
that, in the Kondo regime, the impurity spin forms an effective singlet with all the spins inside the cloud 
||3l and that, only in this regime, one can always choose / so that ^ may be made comparable with A'^; at 
variance, in the dimer regime, the impurity in \ GSi) picks out - no matter what the values of J' and A'^ are- 
only an individual spin in the chain to form a singlet ( a valence bond) while the remaining spins form 
singlets (local dimers) with their nearest neighbors [I j. Thus, only in the Kondo regime, one may use 
the remarkable resource to select - for any A'^- an initial state \GSi) which is free from local excitations: 
in fact, when = N — 2 (f = J'^pj) there is only a 5'ingle entity, namely the A'^th spin, interacting with the 
impurity-cloud composite. As this situation can always be engineered by choosing, for any A'^, J' = J'^^,, 
this allows for an explanation [8] of the distance-independent entanglement in Fig. [5j). The proposed 
scenario provides an intuitive grasp on why, only in the Kondo regime, an optimal entanglement between 
the ending spins may emerge from quench dynamics as the result of the interplay between very few 
states. At variance, in the dimer regime, the energy released by quenching is dispersed over the variety 
of different quantum modes associated to the local dimers and no significant long range entanglement 
may then be engineered. 

To provide a more quantitative analysis, one may expand \ '^f{t)) in terms of eigenvectors of Hp. By 
exact diagonalization (up to A'^ = 20), one finds that, in the Kondo regime, only two eigenstates of Hf 
(the ground state j^i) and one excited state |£'2)) predominantly contribute to the dynamics: 



f,00\ 



\Ex) = ai\\if-)iN\r)h + M^^)iN\r 

+ \n)M<i>'%-\w^W)h) 

Wi) = a2\w-)iN\r)b-HmvN\^'''')b 

+ \n)m\^'')b-W^)\^^)b)- (4) 



In Eq. dll) the first and the last spin are projected onto the singlet {\\lf^)) and the triplets (|00), 1 11) and 
|VA+)) while the states of all spins in the body of the chain have been specified by the index b. After a 
time t the state evolves as 

W{t)) = (£i|GS/)|£i) +e-'^^'(£2|G5/)|£2) + ..., (5) 

where, AS is the energy separation between the two levels. One defines t = topt as the time for which the 
contribution of \ Y^)\N\'P^)b is most enhanced in \ \l/{t)) due to a constructive interference. 
The condition for the onset of constructive interference is 



\{E,\GSi)p,\^\{E2\GSi)l52\, (6) 

so that terms other than \ Y^)iN\<l>^)b in \¥{f)) do not contribute at t = topt- When 7' — )• 1 (very small 
cloud) then KSijGS'/)! 1 while |(£'2|G5'/)| « as the ground state is hardly changed on quench: thus, 
constructive interference between j^i) and \E2) it is impossible. The condition of Eq. Q cannot be 
satisfied also when /' — 0, since one has now that jSi ~ and j32 ~ 1 (the end spins form a singlet and 
triplet with each other in l^i) and \E2) respectively |!2T||). Thus, only for intermediate /' entanglement 
may peak. When J' > Jopt < N — 2)- and particularly for ^ < N/2- the entanglement between the 
ending spins is frustrated by the existence of local excitations whose number increases as the size of 
the cloud gets smaller. In addition, when /' < J apt {£, > N — 2), the Kondo cloud overtakes the chain 
and the A^th spin is already included in the cloud and its tendency is to screen the original impurity as 
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Figure 6: (a) t„pt vs. 1 / -s/J^t the Kondo regime, (b) £"„, at ? oc 1 //^^^^^ vs. J2 for chains of different 
lengths, (c) Entanglement attained at tgpt vs. J' for Kondo and dimer regimes for A'^ = 20. (d) £,„ vs. 
temperature after bond quenching (blue line) and induced by (static) weak coupling with the rest of the 
chain (red line) for a chain of N = 10. 

in \ Y^)\N\'l>^)b rather than to pair with it to form a spin one (as in the last three terms of l^i)). This 
makes /3i quite small, and it becomes smaller as the cloud overtakes the chain and again the condition 
of Eq. ^ cannot be fulfilled. Consequently, the optimal situation is realized when J' = J'^p^ (t, = N — 2 
), i.e just before the cloud overtakes the chain. Thus, only in the Kondo regime, one can convert-for any 
A^- the useless entanglement between the impurity spin and the Kondo cloud into a usable entanglement 
between the ending spins of the chain. The emerging long distance entanglement analyzed in this paper 
is, indeed, a genuine footprint of the presence of the Kondo cloud in \GSi) (Fig. Efa)). 

In Fig. [6J5 we plot- for both regimes- the entanglement reached after waiting for a time interval of 
the order of 1 /J'jpf One notices that, for J2 > entanglement peak decreases sensibly and goes to 

zero rather soon. The plot of the maximal entanglement vs. /' is given in Fig. [6]:: here one sees that, in 
the Kondo regime, the entanglement rises from zero already at very small values of J'. This is expected 
since, in the Kondo regime, to a small J' is associated a large cloud containing the impurity sitting on the 
left side. 

The essential role of the Kondo cloud in the entanglement generation is further probed if one let 
evolve the ground state \GSj) with a doubly quenched Hamiltonian obtained from Q by isolating (i.e., 
putting /' = 0) the left hand side impurity while keeping fixed to J'gpi the bond connected to the last spin 
(see Fig. |7J)). This forbids the dynamical build up of that "portion" of the total entanglement which is 
only due to end to end effects. In Fig. [8] we have plotted the vs. N after double quenching in both 
regimes. Fig. [8]shows that entanglement in the dimer phase collapses already when N > 12 while it stays 
unexpectedly high - and almost independent on A/^ - in the Kondo regime; the existing entanglement of 
the Kondo cloud with the impurity |j3l is dynamically swapped over to the last spin. 

3.3 Thermal Robustness of the Kondo Cloud Mediated Entanglement 

Note that a long distance singlet between the end spins may be realized in a ground state when those 
spins are very weakly coupled (/' = s/y/N « \/^/N) to a spin chain [21]. This static approach to 
generate entanglement relies on couplings which ai^e so weak that they can merely be regai^ded as per- 
turbations. Such entanglement is not robust against thermal fluctuations due to the smallness of the gap 
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Figure 7: (a) Different \GSi) for entanglement generation through quench dynamics. Top: The ground 
states with <N/2 (no entanglement). Middle: N/2 < ^ <N — 2 (some entanglement). Bottom: 
^ = N — 2 (optimal entanglement), (b) Decoupling the impurity from the chain. 




Figure 8: vs. A'^ after decoupling the first impurity. 



(oc 7'2 _ £2^^-) between the ground state and a triplet state between the end spins. On the other hand our 
approach enables to generate entanglement dynamically even for /' as high as 7^^, l/(logA'^)^. Even 
when temperature is increased, the entanglement is not seriously disrupted till K^T exceeds the Kondo 
temperature l/t, = \/{N — 2)) after which the Kondo cloud does not form. As a result, while in 
the dynamical approach KbT < 1/{N — 2), in the static approach one has KbT < e^/N: thus, the long 
distance entanglement generated through the dynamical approach is thermally more stable. For instance, 
for £ ~ 10^ \ our dynamical approach is robust for temperatures 100 times higher than those required 
for the static approach. In Fig. [6]l, we plot £"„,- as obtained in both approaches- vs. temperature for 
A'^ = 10. In the static approach, the ground state is replaced by the thermal state, whereas in our dynamic 
approach it is the initial state which is taken to be the relevant thermal state. We ignore thermalization 
and relaxation during dynamics since the dynamical time scale, set by topt, is fast enough (this is also an 
advantage over slow dynamical schemes with weak couplings 1221 ). 

We should also point out some systems where the Kondo cloud mediated long distance entanglement 
may be observed such as spin chains in ion traps ll32l . with trapped electrons |[33l . in chains of P donors 
in Si ll34l and Josephson chains with impurities (351. 
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4 Summary 

To summarize, we reviewed our analysis |[3l of the ground state entanglement of the Kondo spin chain 
model from the viewpoint of a genuine entanglement measure, namely the negativity. This readily 
showed that the impurity spin is indeed maximally entangled with the Kondo cloud; we provided an 
mdependent method to determine the Kondo screening length together with a characterization of the 
ground state of the Kondo spin chain in the Kondo regime. We proposed an improved DMRG approach 
enabling to account for the entanglement between the impurity and a block of spins located at the other 
side of the chain for different lengths of the block. We defined an Entanglement Healing Length EHL 
and showed that, in the Kondo regime, the EHL L* scales with the impurity coupling /'just as the Kondo 
length does. Finally, our approach showed that, in the Kondo regime, the entanglement scales exponen- 
tially with L/L* and that, in the gapped dimer regime, though it is still possible to define an EHL, the 
impurity-block entanglement is usually smaller and has no characteristic length scale. 

We have also shown that substantial Zong range distance independent entanglement can be engineered 
by a non-perturbative quenching of a single bond in the Kondo regime- and only here!- of a Kondo spin 
chain. This is the first example where a minimal local action on a spin chain dynamically creates long 
range entanglement. In contrast to all known schemes for entanglement between individual spins in 
spin chains, here the entanglement attains a constant value for long chains rather than decaying with 
distance. This behavior arises since, in the Kondo regime, one may always select the initial state so as to 
fulfill the condition ^ = N — 2. As the coupling is non-perturbative in strength (J'^pj l/(logA'^)^), the 
entanglement generation is both fast and thermally robust. The long distance entanglement mediated by 
the cloud and the periodic dynamics, provides a clear signature of the existence of the Kondo cloud in 
a quantum system with impurities. These features are not seen when the Kondo cloud is absent, e.g., in 
the gapped dimer regime and- more remarkably- in other conventional gapless systems. 
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